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^ ' Abstract 

D ■ Bekenstein has shown that violation of Weak Equivalence Principle is 

strongly supressed in his model of charge variation. In this paper, it is 
, shown that nuclear magnetic energy is large enough to produce observable 

' effects in Eotvos experiments. 

1 Introduction 



The variation of fundamental constants has been an important subjeet of re- 
search since Dirac stated the Large Number Hypothesis (LNH) [D [2] ■ In the 
latter times, the interest in that subject has been aroused again since such a 
variation is a common prediction of several "Theories of Everything" (TOEs), 
such as string theories [3J. One possible low-energy limit of these TOEs is 
I Bekenstein's variable charge model [UISIH], since it has all desirable properties 

that such low-energy limit should exhibit. 
' Since Dirac's proposal, many attempts have been made to detect the pro- 

. posed variations, most of them with null results (for reviews, see [21IZ1IH])- An 

' interesting possibility is that a space variation of fundamental constants should 

0^ . produce a violation of the Weak Equivalence Principle [3], a fact that can be 

proved easily using energy conservation lOj. 



^ . The Weak Equivalence Principle (WEP) states that the world line of a body 



immersed in a gravitational field is independent of of its composition and struc- 
ture [SHU] , a generalisation of Galileo's law of Universality of Free Fall: the local 
' acceleration g of a body is independent of its composition and structure. Since 

General Relativity has the Equivalence Principle as one of its consequences, 
testing for its validity is an important form of the search of "new physics" . 

The most sensitive forms of those tests are the Eotvos experiments: testing 
the equality of acceleration for bodies of different of different composition or 
structure [T^ EH H] . Several accurate tests have been carried in the second half 
of the 20th century and up to now O US [H [T71 [H [H] . These tests impose 
strict bounds on parameters describing WEP violations [HI . 

However, in his 2002 paper Bekenstein 5 proved that a violation of WEP 
is highly unlikely in his model. We shall discuss briefly this issue later on, but 
the origin of this statement is a wonderful cancellation of electrostatic sources 
of the V field, leading to a null effect in the lowest order. No such cancellation 
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happens for magnetostatic contribution, but a simple examination of the Solar 
System magnetic energy density suggests that a breakdown of WEP should be 
inobservable. 

In this paper, we discuss the detection of a space variation of a in Beken- 
stein's model, considering the fluctuations of magnetic fields in quantum sys- 
tems. 

The rest of the paper is structured as follows: in section [5] we make a short 
summary of the main results on Bekenstein's model related to our problem; 
section [3] deals with the motion of a composite nonrelativistic body in external 
gravitational plus dilaton fields, to find an expression for its anomalous acceler- 
ation; in section |4] we discuss the magnetostatic energy of matter in a quantum 
system and in section [5] we state our results and conclusions. The appendix is 
devoted to a simple proof of equation p^. 

2 A survey of Bekenstein model 

Bekenstein's proposal [H [S] was to modify Maxwell's electromagnetic theory 
introducing a field e to describe a variation. An unique form of the theory (up 
to a parameter) was found from the following hypotheses: 

1. The theory should reduce to Maxwell's for a constant a. 

2. a variation must be dynamical (i.e. generated by a field ijj — hie). 

3. The dynamics of the field is derived from a variational principle. 

4. The theory must be causal, gauge and time-reversal invariant. 

5. The smallest length in the theory should be the Planck length ip. 

The latter statement should be dropped if the theory is considered a low- 
energy limit of some TOE, since these introduce other fundamental length scales. 
The application of the above hypothesis lead to a unique form of the action 

S — Scm + + S'mat + Sq, (1) 

where 

Som = J e-^'^rf^.^gd^x (2a) 

is the modified Maxwell action [/^^ = (e'^A,,),^ - (e'^A^),^ = e'^F^^] ; 

= ^ / {d.^ fV^gd^x (2b) 

is the Ip field action and S'mat and Sq are the matter and Einstein actions. The 
local value of the electric charge is 

e(2;^) = eoe^("") ^(x^) = e^^^^'^ao (3) 

where eo,Q;o are reference values of these magnitudes. 
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The general equations of motions for these fields are 



□^ = ffl^-^^), (4b) 



{e'^F^'')^^^i7rr, (4a) 
he \dip Stt 

the latter quantity being the rest mass energy density. 

A word of advice is due: in his papers [U [5] Bekenstein uses an ensemble 
of classical particles to represent matter. This is not a good model of matter 
wherever quantum phenomena are important, neither at high energy scales or 
small distances scales, since fermions have a "natural length scale" , namely the 
Compton wavelength of the particle Ac = h/mc. One must be wary of jumping 
to conclusions in these regimes. (See also [21]). 

From the above equations of motion Bekenstein 5J derives several theorems. 



Cancellation theorem For a electrostatic field equation (|4bp can be written 
in the form 

V • (e-2'''£;) = 4ttp (5a) 



da . n,i,E' 



e 



'2 



dtp 47r 



(5b) 



In the source term for ^ the first term cancels almost exactly the second 
and the asymptotic value of tJj is almost exactly suppressed. 

WEP for electric charges The equation of motion of a system of charges in 
an electric field, in the limit of very small velocities, reduces to 

MZ = QE (6) 

and where M and Q are the total mass and charge of the system. Thus, 
there is no WEP violation. 

The above results use the classical point charges model of matter. On the 
other hand, the equation of motion for ijj for a static system of magnetic dipoles 
is 

VV = -We-2^?^ (7) 
An 

and there is no cancellation of sources. From an estimate of the field intensities 
in the Solar System, Bekenstein states that no observable WEP violation can 
be detected in laboratory experiments. This latter result is also based in the 
classical point charges model of matter. 



3 Motion of a composite body in the tl) field 

Let us now find the Lagrangian of a body composed of point-like charges, such 
as an atom or an atomic nucleus. We shall work in the nonrelativistic limit 
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for the charges, but we shall keep for the moment the full expression for the 
electromagnetic field. We shall treat the system as classical and later on quantize 
it in a simple way. The techniques we use are a lightweight version of those used 
in the TiJe/i formalism [531 E] • We assume that there are external dilaton -0 
and Newtonian gravitationa fields acting over the body, but we shall neglect 
the self fields generated. With these approximations, the Lagrangian of the 
system takes the form 



L 



Mtot VP] 



E 



-mavl ~ ma4>NiXa) " Ca^iXa) -Va • A{Xa) 



(8) 



{E^ - B^)dV 
16n ■ 



To eliminate the electromagnetic fields we use the equations of motion ([4al 
together with the Lorentz gauge condition to obtain 



These equations can be solved using retarded potentials 

Jct(irot) 



R c 



R 



-dV', R=\x~x' 



whose slow-motion approximations are 



' R 



B ^ - 

c 



1 f jcf X R 



i?3 



dV. 



(9) 



(10) 



(11) 



In these equations we shall neglect the contribution of since they are 
much smaller than the usual current contribution. After some transformations 
the Lagrangian can be written as 



L = 



'1 



E 

a 

1 

2? 



rriaV^ - ma4>N{Xa) 



„24> Pc{x)pc{x') 



R 



dVdV' 



(12) 



R 



dVdV'. 



where we have replaced sums over pair of particles or currents with integrals. 

For a macroscopic solid body we shall be interested in the motion of the 
center of mass. The separation of this motion is easily achieved with the usual 
substitutions and developing the slowly varying external fields 4>n and ■0: 



■CM 



Eg IT^aXa 
Mtot ' 

Xa = RCM + x'^. 
(t>Nix) ~ 0Ar(ilcM) - 



VcM = E'""' 



(13) 



Va = VcM+v'^, (14) 

o(x''), i;{x) ^ i^iRcu) + O (^x'^y, (15) 
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and besides 



Mtot m ^ Mtot [0] 



dMc^ 
dip 



iPiRcu) + O 



(16) 



Substitution of the above in the Lagrangian leads to 



''cm 



Oat 



(-Rcm)- 



(17) 



■2iP{Rcm)E„ 



The electrostatic contribution cancels with the mass dependence on -0, ac- 
cording to Bekenstein theorem, and the neglected terms are of either tidal order, 
negligible in laboratory tests of WEP, or of higher order in t/i. 

The above Lagrangian shows that a body immersed in external gravitational 
and Bekenstein fields will suffer an acceleration 



Em 

RcM = a = g + 2-^ VV'IcM • 



(18) 



The latter term is the anomalous acceleration generated by the Bekenstein 
field. The acceleration difference (|18p is tested in Eotvos experiments. 



4 Magnetic energy of matter 

In a quantum model of matter, magnetic fields originate in not only in the 
stationary electric currents that charged particle originate and their static mag- 
netic moments but also in quantum fiuctuations of the number density. These 
contributions to the magnetic energy have been computed in Ref. [22l [11] from 
a minimal nuclear shell model. The matrix elements of the current operator can 
be related to the strength of the dipole resonance, with the result 

where R{A) is the nuclear radius, E is the giant dipole mean absorption energy 
and J adE its integrated strength function. These quantities have the following 
approximate representation 

R{A) = 1.2yl3 fm, E - 25 MeV, j adE ~ l.GAMeV fm^. (20) 

Since the magnetic energy density is concentrated near atomic nuclei, it can 
be represented in the form 

eM ^Y^^rJix - Xa) ^Y.^rnnb{x), (21) 

a b 

where index b runs over different nuclear species. Define 



as the fractional contribution of the magnetic energy to rest mass. Then 

e,n{x) = (,n{x)p{x)c^, (23) 
where p{x) is the local mass density and 

Ux) = (24) 

is the local mass-wheighted average of • 

With expression we can write equation ^ in the form 

V^V' = -STTK^c^e-^'^Cmp. (25) 

For small -0 we can find a solution for an arbitrary distribution of sources 

1 r - 

■0 — Sttk^c^— / x'^C,m{x)p{x)dx (26) 



whose asymptotic behaviour can be expressed in terms of the newtonian gravi- 
tational potential 

^ 'GM^^^^' Jo ^ Cn.i^)P{x)dx = 2 ( — 1 Cm— (27) 

where Cm is the mass averaged value of Cm and we have introduced the Planck 
length £p. 

5 Results and conclusion 



From ([27|) we obtain for the differential acceleration of a pair A, B of different 
bodies 

,„.,s, = ^i^ = 4(^£)\.,(a-ce) = c,(^s)' m 

where Cs:Ca7C-B are the magnetic energy fractions of the source, body A and 
body B respectively. 

Ej_ 

Mic^ 

Table [T] shows the results of the most accurate versions of the Eotvos exper- 
iment. A simple least squares fit with the statistical model y — CfX^ yields 



^ , 0.0003 ±0.0006 (30) 

Jp 

from which we get the "3(t" upper bound 

<0.002 -f- < 0.05 (31) 

tp J tp 

This last equation encodes the main result of this paper: strict upper bounds 
can be set from Eotvos experiments on the Bekenstein parameter ^b/^p even if 
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A 


B 


Source 




10^^77(A, B) 


Ref. 


Al 


Au 


Sun 


17.5 


1.0±1.5 


[Til 


Al 


Pt 


Sun 


17.5 


03±0 045 


[15] 


Cu 


W 


Sun 


8.8 


0.0±2.0 


[16] 


Be 


Al 


Earth 


6.8 


-0.02±0.23 


[TT] 


Be 


Cu 


Earth 


10.4 


-0.19±0.25 


[TT] 


Be 


Al 


Sun 


16.1 


0.40±0.98 


[H] 


Be 


Cu 


Sun 


24.6 


-0.51±0.61 


[H] 


Si/Al 


Cu 


Sun 


8.8 


0.51±0.67 


[17] 


EC 


MM 


Sun 


-7.6 


0.001±0.032 




Be 


Ti 


Earth 


6.9 


0.004±0.018 


m 



Table 1: Results of Eotvos experiments. The columns show the composition of 
the bodies, the source, the coefficient of (^s/^p)^ in equation (|28p . the measured 
value of rj and its la error. 

the electrostatic field does not generate tj) field. These bounds are much larger 
than the ones that would result if electrostatic energy density would generate 
■0 field intensity. This calculation was carried in the 1982 paper of Bekenstein 
[3] and has been repeated several times (e.g. [Ml HOI 123) with the result 

(g)^^< 8.7x10- (32) 

one order of magnitude smaller than (|3ip . 

It is interesting to compare our result pip with the results obtained from an 
analysis of all evidence from time variation of the fine structure constant a [25j . 
In that paper, an effective value of C = 10~^ was used, following the suggestion of 
ref. [3] and a la bound on {iB/^pY < 0.003 was found. From the estimate of 
in reference [S] we compute an effective value oi (u = 2.7x IO^^Hb — 1.4x 10^^ 
and so we find a 3a upper bound 

/ < 0.8 (33) 

one order of magnitude larger than pip . 

In conclusion, we have shown that very strict bound can be put on the 
Bekenstein model parameter Ib / I^p from the quantum fluctuations of the mag- 
netic fields of matter. From equation pip one should discard the Bekenstein 
model, but since it can be obtained as a low energy limit of string models, the 
latter conclusion should be taken with a grain of salt. 
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A Proof of ([19]) 

Reference [22] does not give a proof of equation p9p . The following proof is 
based on their methods. 
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Let us write the total magnetic energy of the nucleus in the form 

a 

where a runs over a complete set of eigenstates of the nuclear hamiltonian H . 
The current operator is defined as 



j{x) ^y^5{x- Xa)ea — , (35) 



where the sum runs over all particles in the system. Neglecting the momentum 
dependence of the nuclear potential, we can write 



= - [Xa,H\ 

rria n 



Substitution of the above in equation ([35]) yelds the result 
{0\jix)\a) = ■^^J(a;-a;a)(£;o-Sa)(0|e,a;„| 

a 

= {X - Xa) {Eo - Ea)doa, 



(36) 



with d{x) the polarization (dipole density) operator. 

If we assume a constant density within the nucleus, the dipole density can 
be represented as 

Vn 

where Vat = is the nuclear volume, and so 

(0| j(a.) |a) . {a\j{x) |0) ^ \^^cos9, (37) 

where 9 is the angle between x and x' . Thus, the magnetic energy can be 
expressed approximately as 



^" - — vl, • ^^^^ 



The last factor is equal to The first one can be computed from the 

connection between the strength function and the photoabsorption cross section 

47r 2 

(Joq = 1— £^qoMqo| • (39) 
he 

From this, we easily get 
where _E ~ 25 MeV is the mean absorption energy, roughly independent of A. 
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The cross section satisfies the Thomas-Reiche-Kuhn sum rule 

/ (j{E)dE = (1 + ^ - (1 + x)15 MeVmbarnA, (41) 

J mc A 

where x ~ 0.2 takes into account exchange and velociy dependence of nuclear 
interactions. Combining equations ([55]) . ([3D]) and (HT]) we obtain equations ([TO]) 
and (PUj) . 
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